On the theory of ferroelectric sohd solutions 



A.Yu. Zakharo\{!| and M.I. BichurirQ 

Novgorod State University, Veliky Novgorod, 173003, Russia 
S. Priya 
Virginia Tech, USA 

The paper contains an application of lattice model to ferroelectric solid solutions. Short-range 
parts of interatomic potentials are taken into account by means the lattice structure introduction. 
Long-range parts are considered in effective field approximation. The equations for the mean values 
of the components contributions into polarization are derived. The general equation for the Curie 
temperature dependence on solution composition is obtained. Theoretical calculation of Curie tem- 
perature of the system Pbi_a;Sna;Ti03 fulfilled. 
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I. INTRODUCTION 

Let us consider a ferroelectric solid solution as a lat- 
tice with two kinds of dipoles distributed over the sites. 
Suppose the energy of two dipoles located in sites Rj and 
R^- has the following form 1, 2] 



where 



Urj = - ^ (R. - Rj) (D« (R.) • (R,)) 

a, 7 

xua (Ri) (Rj) , 



(1) 



where Ua (Ri) is a dichotomous random variable with 
values 1 and 0: it is equal to 1, if site Ri contains a 
particle of a-th component, and otherwise. 

These variables ni (Ri) and n2 (Ri) obey the relation 

ni(R,)+n2(R,) = 1, (2) 

their mean values are 

(ni(Ri)) = ci, (n2(Ri)) = C2, ci + C2 = 1 (3) 

where Ci , C2 are components concentrations in the solu- 
tion. 

The functions Q^j (R) in lattice models [1, Q are the 
long-range parts of interaction potentials with cutting 
on short distances. A system of dipoles at presence of 
external electric field E has the following Hamiltonian 

^ = E E (R» - Rj) (Da (RO • (R,)) 



„ (R,) (Rj) - ^ (E • D„ (R,)) n„ (R, 



(4) 



In self-consistent field approximation Jd\ for the long- 
range parts of the interactions this Hamiltonian trans- 
formed to following form: 



(5) 



(Ef • (R,)) Ua i^R^) 



Ef =E 



"•"O ( E '3a7 (R« - Rj") D7 (Rj) (Rj 



\ J, 7 

The last term 



(6) 



(7) 



^ ( E (R« - ) D7 (R, ) (R, ) ) (8) 



in this formula is the local effective field at point Ri due 
to all the dipoles for dipole of a-th kind (in general case 
the effective fields for the components are not identical). 



II. POLARIZATION EVALUATION 

A. Effective fields 

To effective fields (O evaluate it should to find the 
mean values {n^ (Rj)) under the sum over j, 7 sign in 
In general case these mean values depend on the correla- 
tions between components distributions, but in effective 
field approximation Q these mean values are concentra- 
tions of the components: 



(«7 (Rj)) = C7 



(9) 



As a result we have the following expression for effective 
field 



Ef = E+i VQi°^(D,) 



2 ^ 

7 
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(10) 
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where 



(the last relation holds by virtue of translation invariance 
of crystal). 



B. Generating functional and polarization of 
solution 



Generating functional (i.e. the partition function a 
system at an external field presence) in effective field ap- 
proximation has the following form 



1 



X exp 



47r 



^E (Ef -D^ (R,))rv (R, 



(12) 



where dVli^a is a solid angle infinitesimal element in spher- 
ical coordinates {Da, 9i,^i) 



sin 6i dOi d(pi , 



(13) 



/3 — l/T, T is absolute temperature in energetic units, 
iVi and N2 are numbers of dipoles. 

The multiple integral ([T^ is a product of the same 
type integrals 



y - exp [/.i^-i. cos 



(14) 



therefore generating functional in mean field approxima- 
tion can be evaluated easily 



Z 



1 



Ni\N2l 



sinh {(3EfDi 



Ni 



sinh {(3EfD2) 



A'2 



(15) 

Hence we find the system polarization (i.e dipole moment 
per volume unit) 



P 



1 d\nZ 



V 



N2 
' V 



VP dE 



coth {PEfDi 



1 



coth {PEfD2) - 



{PEfDi 
1 



(16) 



{PEfD2 



Do 



Thus, for binary solid solution full polarization con- 
tains two contributions related two of the components, 
but these contributions are mutually dependent via ef- 
fective fields. Therefore the components contributions 
into polarization are not additive. 



C. Mean value of a single particle dipole moment 
in effective field 



Mean value of a single particle dipole moment in ex- 
ternal field E defined as: 



{D) = 



J e-'3^i D cos6'df}i 



and has well known result 

(D) = DL{f3DE), 

where 



L{z) = 



coth (z) 



(17) 



(18) 



(19) 



is the Langevin function. Hi is the Hamiltonian of single 
dipole with moment D in external field E 



77i = - (E • D) . 



(20) 



For one-component system expression ([7]) has a simple 
form 



E<=ff = E i Q(") (D) , 



(21) 



Substitution E°^ instead of E into formula (IT51) leads 
to transcendent equation with respect to {D). It is well 
known this equation has nontrivial solution {D) ^ at 
T <T,. 



D. Polarization of solid solution 

System of equations for mean values (-Di) and {D2) of 
dipole moments of the solution components follows from 
relation (IT6l): 



{Di)=DiL\(iDAE 



{D2) =D2L{I3D2{E 



7/ ^^7 



7 



(22) 

Note the interdependence between (Di) and {D2) realizes 

via non-diagonal elements of matrix Qaj only. 

It is clear that the solution {{Di) , {D2)} of this sys- 
tem of equations depends on temperature, external field, 
and all the matrix elements Qa-y- To solve this system 
we should know the parameters of the physical system, 
i.e. matrix elements. This elements should be find using 
some experimental data. 

Solution of system ((22|) with respect to (Di), {D2) per- 
mits to find polarization of the system as function of ex- 
ternal field E 



N 



(23) 
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and consequently susceptibility x of this system: 



iV 



■ d{D,) 



d{D2) 



(24) 



The Langevin function L{z) contains two terms. The 
first term is algebraic, the second is the transcendent 
term. Both of them have a singularity at z = 0. These 
circumstances complicate search of the solution. There- 
fore, the Langevin function should be approximated by 
some more suitable function with correct asymptotic be- 
havior at z — ?> and z — > oo. In the capacity of such 
approximation we shall use the following function 



L{z) w — arctan f — ) 

TT V 6 / 



(25) 



Note, this approximation is well not only for Langevin 
function, but also for its derivative 



L'iz) 



12 



36 -I- 7r2z2 • 



(26) 



Differentiating both of sides of equation ([28)1 with re- 
spect to E 



dE 



12/3 [D^f 



X ( i + in(o)^i^ 



(30) 



we find derivate 



12/3 (Do 



dE 



- 6/3Qi°2 {D^) 



(31) 

Eliminating Q^aa from the last equation with ac- 
count ([29]) . we obtain susceptibility of the pure ferro- 
electrics: 



X ■ 



d{D^) 
V dE 



This approximation for Langevin function permits to 
simplify the system of equations (p2)) 



2 Da, ttPD^ 
arctan — - — 



E 



2 ^ 

7 



(27) 



V 



12/3 (D„)' 



36(l-^) + (f)' 



ED + 3T, 



{-De 



-De 



(32) 

At -E = the susceptibility have a singularity in vicinity 
of the critical point Tc- 



III. PARAMETERS OF THE COMPONENTS 
A. Curie temperatures of the components 

Some of the parameters Da and of solid solution 
can be find from experimental data of the components. 
Equation (l?fl) for a pure component (cq = 1, c-y\^_^^ — 0) 
have the following form 



2 Da 



X arctan 



TTfiDg 
6 



(28) 



This equation describes connection between external 
field E and mean value of dipole moment Da- Critical 
temperature can find from condition existence of non- 
trivial solution (i.e. (Da) 0) in external field vanishing 
E — 0. Graphical analysis of the equation ([TS]) leads 
to following connection between critical temperature 
and model parameters 



{Daf Q 



(0) 



(29) 



B. Polarization of pure components 



Let us introduce the new variables in equation 



Da 



EDg 

T 

T 
TV' 



, < 1; 



(33) 



These dimensionless variables have the following physical 
interpretation: Xa is ratio mean value of dipole moment 
(Da) to absolute dipole moment Da', s (dimensionless 
external field) is ratio the dipole energy E Da in external 
electric field E to temperature T; r is ratio the tempera- 
ture T to the Curie temperature Tc defined by ([29]) (thus, 
T is dimensionless temperature). 

Equation (|28p in these variables determines connection 
between e, t, and Xa- 



6, /TT \ Xa 

— tan [ — Xa ] ~ o — . 

TT V2 / T 



(34) 



Connection between dimensionless external field e and 
mean value of dipole moment Xa at fixed values of tem- 
perature r are presented on Figs.[T] — [SJ 
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FIG. 1. Calculated dependence of the dipole moment mean 
value Xa on external field e at r = 0.3. 



FIG. 3. Calculated dependence of the dipole moment mean 
value Xa on external field e at r = 0.9. 




FIG. 2. Calculated dependence of the dipole moment mean 
value Xa on external field e at r = 0.6. 



At T > Tc (t > 1) the curves Xa (e) pass through the 
origin. Otherwise (i.e. at t < 1)) there are some sponta- 
neous polarization due to nonzero values of equilibrium 
mean dipole moments at £ = 0. Dependence of dipole 
moment equilibrium value Xa on temperature r follows 
from equation (|34p at e = 0: 



CKXaX TTXa 



0. 



(35) 



Dependence of the spontaneous value of mean dipole 
moment Xa on temperature r presented on Fig. [51 
Spontaneous polarization P of pure ferroelectrics ex- 
presses via Xa'- 



P 



N 
V 



(36) 
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FIG. 4. Calculated dependence of the dipole moment mean 
value Xa on external field e at r = 1.2. 



where y is number of the dipoles N per unity of the 



system volume V . 



IV. SUSCEPTIBILITY OF SOLID SOLUTIONS 



After differentiating both sides of each equation in ([T7)) 
over E we obtain a system of linear algebraic equations 
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This system of equations has the fohowing short form: 



dE 



dE 



where 
Aa - 

and 



12/3 



2^a7^7- 



(38) 



(39) 



(40) 



FIG. 5. Calculated dependence of the dipole moment mean 
value Xa on external field e at r = 1.5. 
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FIG. 6. Mean value of the dipole moment function on 

dimensionless temperature t < 1. 



Hence we have for ^^-^ and 



dE 



Ai (1 - A2B22 + A2B12) 



1 — AiBii — A2B22 + A1A2 [B11B22 — B12B21] ' 



d{D2 
dE 



A2il-AiBn+AiB2i) 



1 — AiBii — A2B22 + A1A2 [B11B22 — B12B21] 

(41) 

This system of equations will be used for critical point 
finding. 



V. CRITICAL POINTS OF SOLID SOLUTIONS 
AND MODEL PARAMETERS 



The solution susceptibility (|24p in the critical point has 
a singularity due to vanishing denominators in right hand 
sides of dH]) 

l-AiBn-A2B22+AiA2[BiiB22 - B12B21] =0. (42) 
This equation takes place under the conditions 



for 



d{Da] 

dE 



^ dE 



36+(^/?Z?„)' Qi°^c^(D^) 



(37) 



= 0; 

{D,) = (D,) = 0. 



(43) 



The second of these conditions due to polarization van- 
ishing at T > Tc. 

Substituting expressions (15^ and (HH)) with condi- 
tions account into equation (jH]), we obtain the 
quadratic equation with respect to critical tempera- 
ture T: 



^2 {DifQ^'^c, + {D2rQ)2^c2 ^ 



'^-§-{D,f{D2f 



Vll V22 ^12 V21 



(44) 
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The discriminant of this equation in relation to sym- 



metry property of the matrix Q^a'y is non-negative 



+4C1C2 {D.fiD^fQ^^^Q^^^ > 



(45) 



therefore this equation has two real solutions. 

Two variants are possible depending on the quantity 
Q sign 



(46) 



1. Q < 0. In this case the only of the solutions is 
positive, and the second solution is negative. The 
positive solution has a physical sense as Curie tem- 
perature, the negative solution has not any physical 
sense. 

2. Q > 0. In this case both of the solutions are pos- 
itive. A physical interpretation of this case can 
be clarified after temperature analysis of the equa- 
tions (P7)) solutions. Thus, at condition Q > 
two-component ferroelectric solid solution has two 
critical points. 

Note that the critical points of solutions depends on 
such combinations dipole moments Da and interatomic 

potentials Qa] 

Ga^^Ql^^DaD^. (47) 
Transform the equation (j44p using these combinations 
G22X 



y2 Gil (1 - ■^) 



■T- 



x{l — x) 
36 



G = 0, 



where 



(48) 



(49) 



X = C2 is second component concentration, ci = 1 — a;. 

The left hand side of the equation (pS)) is linear func- 
tion with respect to Gn, G22, G, therefore the simplest 
way of the parameters Gn, G22, G finding from experi- 
mental data is the least squares method. 

In order to realize this method, let us introduce the 
function / (G, Gn, G22, T, x): 



/(G,Gii,G22,T, x) 



Gxx{\ - x) + G22X ^ , x{\-x) 



(50) 
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36 



-G. 



Then the parameters G, Gn, G22 can be find by function 
-F (G, Gil, G22) minimization: 



F (G, Gil, G22) 
{G,Gii,G22,Ti,Xi) = min. 



(51) 



where Ti the experimental points. 

Minimization of this function leads to a system of lin- 
ear algebraic equation with respect to unknown parame- 
ters Gil, G22, G and this system can be solved easily. 

There is the second way. As a rule, the accuracies of 
the experimental data for Curie temperatures at points 
X — and a: = 1 are higher as at intermediate val- 
ues of concentrations x, hence values of the parameters 
Gil, G22 should be find by the experimental data of 
the Curie temperatures of the components using equa- 
tion (I2H) (with relation (|T7)) account). The last of the 
model parameters G can be find using few experimental 
data of the Curie temperatures Tc of solid solution at 
intermediate values of the concentrations. 

As an example, let us consider the solid solution 

Pbi — X 

SuajTiOa. 



VI. SYSTEM Pbi-:,Sn,Ti03 

There are the experimental research of the system 
Pbi_2;Sna;Ti03 in Virginia University. The measure- 
ments results presented partially in the table 



X 


0.00 


0.50 


1.00 


T(x) 


393 K 


623 K 


763 K 



The first line contains concentration of SnTiOa in the 
system, the second line contains the corresponding Curie 
temperatures in the system. 

Unfortunately, the numerical values of the Curie tem- 
peratures found from graphics have essential inaccuracies 
at least order of 10 K. In addition, there is most accurate 
value of Curie temperature for pure PbTiOa. 

Results of the parameters calculations are: 



Gil = 2358 K 
G22 = 4578 K; 
G = -4.043 -10^ 



(52) 



Results of these parameters using to Curie tempera- 
tures of the solutions are presented on Fig. [T) These re- 
sults are in close agreement with available experimental 
data. 



VII. CONCLUSION 

This paper contains the following results. 

1. Deduction of the closed system of equations with 
respect to mean values of the elementary dipole 
moments of the components in ferroelectric solid 
solutions in effective field approximation for long- 
range parts of inter-dipole potentials. 

2. The analysis of the connection between polariza- 
tion, temperature, and external field strength of 
pure components is fulfilled. These results are pre- 
sented in analytic and graphical forms. 
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FIG. 7. Calculated dependence of the Curie temperature (°C) 
on the composition x of the system Pbi-a^Sn^^TiOs. 
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